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OJTc ABSTRACT 

^3 ' We examined effects of the ultraviolet background radiation (UVB) on the colour- 

. magnitude relation (CMR) of elliptical galaxies in clusters of galaxies in the hierar- 

chical clustering scenario by using a semi-analytic model of galaxy formation. In our 
model the UVB photoionizes gas in dark haloes and suppresses the cooling of the 
■ diffuse hot gas onto galaxy discs. By using a semi-analytic model without the effect of 

the UVB, Kauffmann & Chariot found that the CMR can be reproduced by strong su- 
pernova heating because such supernova feedback suppresses the chemical enrichment 
f--*) ■ in galaxies especially for small galaxies. We find that the CMR also becomes bluer 

| because of the UVB, in a different way from the effect of supernova feedback. While 

, the supernova feedback suppresses the chemical enrichment by a similar mechanism 

OO ' to galactic wind, the UVB suppresses the cooling of the hot gas. This fact induces the 

| suppression of the metallicity of the intracluster medium (ICM). In our model we find 

i— I . that the existence of the UVB can plausibly account for an observed ICM metallicity 

that is equal to nearly 0.3 times the solar value, and that in this case we can reproduce 
^ (— I the CMR and the metallicity of the ICM simultaneously. 

^p" 1 , Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 

O ' formation - large-scale structure of the Universe. 
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ical models based on the traditional scenario. This fact that 
the CMR is a metallicity sequence was also confirmed by in- 
vestigating photometric properties of the stellar populations 
of ellipticals (Ferreras, Chariot & Silk 1999). 

In contrast to the collapse/wind model, recent devel- 
opments on both theory and observation of the cosmolog- 
ical structure formation are revealing that objects such as 
galaxies and clusters of galaxies are formed through hierar- 
chical clustering of smaller objects. Kauffmann & Chariot 
(1998; hereafter KC) applied their semi-analytic model to 
the problem of the CMR and reproduced the observed CMR 
when the model includes the chemical evolution process and 
strong feedback to interstellar media by supernovae. More- 
over, they found a tight luminosity-metallicity relation and 
reproduced other properties of cluster ellipticals such as the 
line indices and the Faber- Jackson relation. 

KC interpreted why the CMR is reproduced when the 
feedback process is strong as follows; they considered that 
giant ellipticals must be formed by mergers of larger spi- 
ral galaxies and that dwarf ellipticals must be formed from 
smaller spirals. Because the feedback strength depends on 
the mass of galaxies, the metallicity of stars in more massive 
progenitors becomes higher than that in smaller progenitors. 



1 INTRODUCTION 



It is well known that elliptical galaxies in clusters of galax- 
^ ies have a tight correlation between their colours and 
magnitudes (so-called 'colour-magnitude relation'; hereafter 
CMR). For example, the rms scatter about the mean CMR 
is typically ~ 0.04 mag in the Virgo and Coma clusters 
of galaxies. It is equivalent to observational errors (Bower, 
Lucey & Ellis 1992). Because it is believed that this tight 
relation reflects formation and evolution processes of ellip- 
tical galaxies, many people have studied this relationship in 
order to understand galaxy formation processes. 

The traditional scenario of formation of ellipticals is 
that a monolithic protogalactic cloud collapses and then 
forms stars over a short time-scale until the formation of a 
galactic wind (Larson 1974; Arimoto & Yoshii 1986, 1987). 
In this framework, Kodama & Arimoto (1997) found that 
the CMR is corresponding to the sequence of the mean stel- 
lar metallicities (the metallicity sequence) by comparing the 
CMRs of high redshift cluster galaxies with those of theoret- 
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Therefore giant ellipticals become redder than dwarf ellipti- 
cals. 

In this work, we introduce the photoionization process 
of galactic gas by the UV background radiation (UVB) in 
our model. We have already investigated the effect of the 
UVB on the luminosity function and the colour distribution 
of galaxies (Nagashima, Gouda & Sugiura 1999), and here 
investigate this effect on the CMR, which has not been con- 
sidered. Moreover, it is also shown how the metal abundance 
in the intracluster medium (ICM) is affected by the UVB. 
Through this work, the difference between the mechanisms 
and effects of the supernova feedback and the UVB in the 
galaxy formation process will be clarified. 

In Section ^, we describe the semi-analytic model used 
here briefly. In Section |ij we show the CMR and show that 
the effect of the UVB seems similar to the effect of the su- 
pernova feedback. In Section ^, we investigate the property 
of the effect of the UVB on the CMR. We will show that 
the CMR hardly depends on the intensity of the UVB if it 
is sufficiently strong, and interpret this result physically by 
using a simple model. In Section 5, the metal abundance of 
the ICM is discussed. Section H is devoted to conclusions. 



2 MODEL 

Here we describe the semi-analytic model that we use briefly. 
At first, merging histories of dark haloes are realized by 
the extension of the Press-Schechter formalism (Press & 
Schechter 1974). We adopt a method given by Somerville 
& Kolatt (1999). The cosmological model is fixed to a cos- 
mological constant-dominated flat universe, Qo = 0.3, J1a = 
0.7, h = 0.7 and erg = 1, where h is the Hubble parameter, 
h = f/o/100km s Mpc , and as is the normalization of 
the power spectrum of density fluctuation. We use a power 
spectrum given by Bardeen et al. (1986). Haloes with cir- 
cular velocity smaller than 40km s _1 are treated as diffuse 
matter. In this paper, we consider haloes only with circular 
velocity V c = 10 3 km s" 1 at z = 0, which is a typical value 
for a cluster of galaxies. 

Next, in the merging path, we calculate the evolution 
of the baryonic component from higher redshift to present. 
Diffuse gas in a newly collapsing dark halo is heated up to 
the virial temperature of the dark halo by shock heating (the 
hot gas). The shock heating occurs only in the case that the 
mass ratio of the largest progenitor among all progenitors 
of the halo to the newly collapsing halo exceeds / re hcat- In 
the other case, the hot gas in the new halo conserves the 
temperature of the largest progenitor. If the new halo has 
no progenitor halo, the mass fraction of the hot diffuse gas 
is equal to flb/fio- Here we use fib = 0.015/i -2 . Assuming 
the singular isothermal density distribution of the hot gas, 
we calculate the cooled gas mass by equating the cooling 
time-scale T coo i(r) to the elapsed time from the last shock 
heating telapse, r coo i(r coo i) = t c ia pS c, where the cooling func- 
tion given by Sutherland & Dopita (1993) is used. If the 
UVB does not exist, the gas within the cooling radius, r coo i, 
cools and is accreted by a disc of the halo central galaxy. 
In the other case, the UV photons penetrate the cooled gas 
from r coo i to ruv, where r coo i > ruv, and the outer layer 
r > ruv is photoionized. This ruv is estimated by assuming 
the inverse Stromgren sphere approximation (Nagashima et 



al. 1999). Thus the gas within ruv is cooled. In this pa- 
per, we assume that the intensity of the UVB, J, evolves 
as J tx (1 + z) 7 , where 7 = 4 for z < 2 and 7 = —1 for 
2 < z < 6. The UVB does not exist before z = 6. We found 
that the results in this work does not depend on the shape 
of the intensity evolution of the UVB qualitatively. In the 
following, we parameterize the UV intensity at z — 2 by 
J_2i = (J/10" 21 erg cm" 2 s" 1 Hz" 1 sr" 1 ). It should be noted 
that in the case where the circular velocity of haloes exceeds 
500 km s _1 , we prevent the cooling process manually. This 
is the same procedure as Kauffmann, White & Guiderdoni 
(1993). 

Stars are formed from this cold gas. The star formation 
rate in discs is given by M* = M co i d /r*, where M* and M co i d 
are the masses of the stars and cold gas, respectively, and r* 
is the star formation time-scale. We adopt a similar form to 
that given by KC, r* = r» [r dy n/Tdyn(z = 0)], where r dyll is 
the dynamical time-scale of the halo which is calculated by 
the spherical collapse approximation (Tomita 1969; Gunn & 
Gott 1972), and r° is a free parameter, which is set to 2Gyr 
in this paper. 

With star formation, supernovae occur and heat up the 
surrounding cold gas to the hot phase (supernova feedback). 
The reheating rate is given by M re heat = f3(V c )M*, where 
P(Vc) = (V c /U h ot)~ ahot , Uhot and a hot are free parame- 
ters, respectively. This is the same parameterization given 
by Cole et al. (1994), and we adopt ahot = 2 according to 
KC. 

Chemical enrichment is solved in a consistent way to the 
above star formation process. The yield y is assumed to be 
twice the solar value according to the strong feedback model 
of KC. While KC assumed that a fraction / of metals is 
released from ejecta of the supernovae to the hot gas directly, 
we do not assume such a process. The metals are released 
to the hot gas by the supernova feedback in the same way 
as for the cold gas. 

We recognize a system consisting of the stars and cooled 
gas as a galaxy. Hereafter we define the mass of a galaxy as 
the mass of stars and cold gas of the galaxy. When two or 
more dark haloes merge together, there is a possibility that 
galaxies contained in progenitor haloes merge together. We 
define a central galaxy as the central galaxy of the largest 
progenitor halo. Other galaxies are denned as satellite galax- 
ies. When the elapsed time of a galaxy from the epoch at 
which the galaxy becomes satellite exceeds the dynamical 
friction time-scale of the halo, the galaxy merges with the 
central galaxy of the halo. When galaxies merge together, if 
the mass ratio of the smaller galaxy to the larger galaxy ex- 
ceeds /bulge, starburst occurs and all cold gas is consumed. 
The same feedback law as that for discs is adopted. Then, all 
of stars becomes bulge stars. In the other case, the smaller 
galaxy merges into the disc component of the larger galaxy 
without any additional star formation activity. In this paper 
we adopt /bulge = 0.2. 

Finally, we calculate the colour and luminosity of each 
galaxy from star formation history of each galaxy by us- 
ing the stellar population synthesis technique. We use a 
simple stellar population model given by Kodama & Ari- 
moto (1997). Through the above procedures, we obtain a 
colour-magnitude diagram of galaxies. In the follows, we pick 
out galaxies with the _B-band bulge-to-disc luminosity ratio 
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larger than 1.52 as elliptical galaxies according to Simien & 
de Vaucouleurs (1986). 

The above procedure is mainly based on Kauffmann et 
al. (1993), Cole et al. (1994, 2000) and Somerville & Primack 
(1999). Though the details of our model are slightly different 
from the previous ones, we confirmed that the main results 
of the previous work, such as both of the field and cluster 
luminosity functions, colour distribution and cold gas mass 
fraction of galaxies, are almost reproduced by our model 
within similar range of values of the parameters (Nagashima 
& Gouda, in preparation). 



3 COLOUR-MAGNITUDE RELATION 

In Figure [lj we show the colour-magnitude diagram for four 
models. The thick solid curves with errorbars indicate the 
CMRs given by the models averaged over 50 realizations. 
The errorbars denote the mean values of la scatter for each 
realization, so they are corresponding to the scatter for each 
cluster of galaxies. Dots indicate ellipticals only for five re- 
alizations. The thin dashed lines indicate the observational 
CMR in the Coma cluster (Bower, Lucey & Ellis 1992) and 
the thin solid lines denote the aperture-corrected CMR in 
the Coma cluster (Kodama, Arimoto, Barger & Aragon- 
Salamanca 1998). Since colour of each elliptical galaxy in 
our model is integrated over the whole region of the galaxy, 
the model CMR should be compared with the aperture- 
corrected CMR. Upper panels show the models without 
the UVB and lower panels show those with the UVB of 
J-21 = 0.1. 

As shown by KC, when we do not consider the effect of 
the UVB, the model with Vhot = 280 km s^ 1 can reproduce 
the observed CMR (Figure gja). Imposing stronger feedback 
with Vhot = 350 km s~ , the colours of ellipticals become 
too blue to be consistent with the observation (Figure |i|d) 
because the metallicity of ellipticals decreases (see next sec- 
tion). The same effect of the bluing is also caused by the 
UVB. In Figure ^p, we show the model with the same value 
of Vhot as that of the model in Figure [j]a and with the UVB 
of J-21 = 0.1. By weakening the feedback strength from 
Vhot = 280 to 220km s _1 , we obtain a CMR in agreement 
with the observation again (Figure |l|d) . This manipulation 
is corresponding to reverse of the procedure from Figure [j^ 
to Qb. Thus we obtain a similar CMR to the observation 
with both of the weaker feedback and the UVB. 

It should be noted that KC found that the CMR reflects 
the metallicity-luminosity relation. We also find that models 
with similar CMR show similar metallicity-luminosity rela- 
tion independent of the existence of the UVB, e.g., the dif- 
ferences of the metallicity-luminosity relations between the 
models of Figures la and Id and between Figures lb and 
lc are negligible. Thus the above properties of the CMR, 
i.e., the dependences on the supernova feedback and on the 
UVB, are related with the stellar metallicity of ellipticals. 
In this connection, we find that the shape of the metallicity- 
luminosity relation of the four models is similar to the metal- 
licity sequence of Kodama & Arimoto (1997), in spite of the 
significant difference between our model and their mono- 
lithic cloud collapse model. This also suggests the impor- 
tance of investigating the chemical enrichment process of 
galaxies. 
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Figure 1. V — K Colour-magnitude diagram. The thick solid 
curves with errorbars denote the model CMRs averaged over 50 
clusters and the errorbars denote the mean values of the la scatter 
for each realization. Each dot indicates an elliptical galaxy only 
for five realizations. The thin dashed line in each panel indicates 
the CMR in the Coma cluster (Bower et al. 1992). The thin solid 
line denotes the aperture-corrected CMR in the Coma (Kodama 
et al. 1998). Model parameters for the four models are shown 
in Table hi Upper panels: J_2i = 0. (a) Vhot = 280 km s — V 
(b) Vhot = 350 km s _1 . Lower panels: J_2i = 0.1. (c) Vhot = 
280 km s _1 . (d) Vhot = 220 km s _1 . Note that the simulations 
presented in (a) and (c) have the same strength of the supernova 
feedback. 



4 CHEMICAL ENRICHMENT OF GALAXIES 

While the above result seems to show that the effect of the 
UVB is the same as the supernova feedback apparently, the 
ways to make ellipticals metal-poor are different. In Fig- 
ure |^a, we show four CMRs with different strengths of the 
UVB, J_2i = 0.01, 0.1, 1 and 10 and with the same feedback 
strength, V ho t = 280 km s _1 . The CMRs are hardly changed 
by the UVB if J_2i > 0.1. On the other hand, the luminos- 
ity functions are clearly affected by the UVB. In Figure ^js, 
we show the cluster luminosity functions of these models. 
The thin solid line denotes the observational cluster lumi- 
nosity function in the Virgo cluster given by Sandage et al. 
(1985). These figures show that the sufficiently strong UVB 
decreases the number of galaxies by its photoionization ef- 
fect while it does not affect the metallicity of each elliptical 
galaxy. 

In the following, we explore the reason why the metal- 
licity does not depend on J_2i- First of all, we evaluate the 
effect of the supernova feedback. The mass-weighted mean 
metallicity of each galaxy is given by (Nagashima & Gouda, 
in preparation) 



(Z*(t))M, = 



^ ^ l-exp[-(l+/3-(fl- Q y)(l-/))^] 
l-exp[-(l + /3-Ji(l-/))^] ' 
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Figure 2. (a) Colour-magnitude relations. Model parameters of 
these four models are the same as those of Figure hjc (Vhot = 280 
km s" 1 ) but for the intensity of the UVB. The solid, dashed, 
dot-dashed and dotted lines denote the models with J-21 = 
0.01,0.1,1 and 10, respectively, (b) Cluster luminosity functions. 
The types of lines are the same as (a). The thin line is the lu- 
minosity function of the Virgo cluster galaxies by Sandage et al. 
(1985). 



where 



F 



7° 

^cool 



1 + 



l + f3-R(l-f) 



(2) 



a is a locked-up mass fraction, a = 1 — R (R is the gas 
fraction returned by evolved stars), and Z® ool is an initial 
metallicity of cold gas at t = t a , where t a is an initial time. In 
this paper, / = and R = 0.25. The starburst corresponds 
to t/r„ — > 00, so the final stellar metallicity depends only on 
the strength of the supernova feedback, /3, as 

(1 - Z° cool )ay 



(Z,(t — > OO)) M, 



1 + (3 - R + ay 



■y _|_ n 1 "cool • (3) 

Therefore the stronger the feedback is, the lower the metal- 
licity is. Thus the differences between Figures |l]a and|l|b and 
between Figures |l|c and |l]d are explained by the strength of 
the supernova feedback, j3. 

Next we investigate the process by which the UVB sup- 
presses the metallicity. When the UVB exists, the amount 
of the cooled gas decreases. Then the amount of the metal 
polluted reheated gas also decreases and the metallicity of 
the hot gas is suppressed. The hot gas, including metal pol- 
luted reheated gas, cools at the later stage again. Therefore 
the metallicity of the re-cooled gas decreases. This corre- 
sponds to decreasing Z® 00 \ at the later stage in eq.(3). Thus 
the mean stellar metallicity is suppressed by the UVB. If the 
UVB is sufficiently strong (J_2i i5 0.1), then Z® ool becomes 
negligible compared with the first term in eq.(3), y/(l + [3). 
As a result, the metallicity does not depend on the J-21. 



5 METAL ABUNDANCE OF INTRACLUSTER 
MEDIUM 

As considered in the previous section, the UVB may affect 
the metallicity of the hot gas. Here we investigate it by a 
simple model. Now consider two ellipticals with the same 
stellar mass which are illustrated in Figure ^. One is exposed 
to the UVB (right in Figure and the other is not (left). 




without UVB 



with UVB 



Figure 3. Schematic descriptions of ellipticals in dark haloes. 
These two galaxies have the same masses of dark matter and 
stars. r v ; r denotes the virial radius of the haloes. The mass of the 
reheated, metal-polluted hot gas is f3M*. The right model has 
weaker supernova feedback than the left model because of the 
UVB. 



Table 1. Metallicity of intracluster medium 



Model 


Vhot (km s 1 ) 


J-21 




Figure la 


280 





0.59± 0.01 


Figure lb 


350 





0.49± 0.01 


Figure lc 


280 


0.1 


0.34± 0.02 


Figure Id 


220 


0.1 


0.38± 0.02 



The second case is required to have stronger feedback (^2) 
than the first (f3i) because the same mass of stars must be 
formed from larger amount of cold gas than in the first case. 
Thus the second elliptical has larger amount of the metal 
polluted reheated gas than the first one. 

From the above discussion, we expect that the metal 
abundance of the ICM of the first elliptical should become 
smaller than that of the second elliptical. We show the 
metallicity of the ICM, Zicm, in Table [j] for the four mod- 
els in Figure |l| Zicm decreases significantly when the UVB 
exists. Because of Zicm — O.3Z0 in some observations (e.g., 
Tamura et al. 2001), the existence of the UVB can plausibly 
account for the metallicity of the ICM. It should be noted 
that KC also obtained Zicm — 0.3 with a non-zero value of 
/ without the effect of the UVB. Their model is nearly cor- 
responding to our model with / rc h oa t — > 0, in which we find 
Zicm — 0.3 with no UVB. Thus our model is not inconsis- 
tent with KC. At present there is no method of determining 
the value of the parameter / rc hcat, although Salvador-Sole, 
Solanes & Manrique (1998) suggest / rc h C at ~ 0.6 by com- 
paring the so-called universal density profile of dark haloes 
(Navarro, Frenk & White 1996) with the theoretical pre- 
diction using both of the spherical collapse model and the 
formation redshift distribution of haloes. Thus we consider 
that models with / rc hcat — 0.5 are more realistic. 



6 CONCLUSIONS 

We investigate the effect of the UVB on the CMR in the hier- 
archical clustering scenario by using a semi-analytic galaxy 
formation model, in which we introduce the photoionization 
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effect of the gas by the UVB. We show that the UVB plays 
a similar role on the CMR to the supernova feedback appar- 
ently. Both effects suppress the mean stellar metallicity of 
each elliptical galaxy. Therefore the CMR becomes bluer in 
the both cases. 

We find that the CMR is affected even if the UV inten- 
sity is weak, J-21 ~ 0.1 and hardly depends on the intensity 
in the case of sufficiently strong UVB, J-21 0.1, in con- 
trast to the supernova feedback. This is because the UVB 
suppresses only the cooling process (see below). Moreover 
we show that the two mechanisms which degenerate in the 
CMR can be distinguished by investigating the metallicity 
of the ICM, and that the existence of the UVB is favoured 
to account for the observed metal abundance of the ICM. It 
will be possible to detect the difference in the metallicity by 
observational studies with XMM-Newton. 

The physical mechanisms affecting the CMR of the 
UVB and of the supernova feedback are different, as fol- 
lows. In the case of the strong feedback without the UVB, 
while much cold gas is reheated before metal-rich stars are 
formed, much metal-polluted gas is expelled in the form of 
hot gas such as the galactic wind. On the other hand, the 
UVB suppresses gas cooling. So the material forming stars, 
that is, the cooled gas, is less than in the former no UVB 
case. Under the UVB scenario, if the supernova feedback 
is weaker than the former no UVB case, a similar mass of 
stars can be formed from a smaller amount of the cooled 
gas, and then the stars are chemically enriched. Note that 
in the case of a UVB the hot gas becomes metal-poor be- 
cause a smaller fraction of metals is expelled from galaxies 
than in the former case. Thus conserving the stellar metal- 
licity of each galaxy by introducing the UVB and weakening 
the feedback, the metal abundance in the ICM can become 
lower compared to the former case of no UVB. 

Recently an overcooling problem has been discussed by 
using high-resolution hydrodynamical simulations (Balogh 
et al. 2001), in which too much gas cools compared with ob- 
servations. The effect of the UVB, as well as strong super- 
nova feedback, will be useful to help solve this problem. Be- 
sides by using the proximity effect, it is suggested J-21 <; 1 
(Jennifer et al. 2000) . In this paper, we showed that the UVB 
affects the CMR and the ICM even in the case of J-21 il 0.1, 
but considered the effect of the UVB only by introducing 
the simple inverse Stromgren sphere approximation. Hence 
these results mean that we should investigate the effect of 
the UVB in more detail. While we believe that we can under- 
stand the effect of the UVB qualitatively, detailed radiative 
transfer calculation will be required for quantitative estima- 
tion of the effect. 
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